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We develop a manifestly Lorentz covariant chiral quark model for the study of baryons as bound 
states of constituent quarks dressed by a cloud of pseudoscalar mesons. The approach is based on a 
non-linear chirally symmetric Lagrangian, which involves effective degrees of freedom - constituent 
quarks and the chiral (pseudoscalar meson) fields. In a first step, this Lagrangian can be used to 
perform a dressing of the constituent quarks by a cloud of light pseudoscalar mesons and other heavy 
states using the calculational technique of infrared dimensional regularization of loop diagrams. We 
calculate the dressed transition operators with a proper chiral expansion which are relevant for the 
interaction of quarks with external fields in the presence of a virtual meson cloud. In a second 
step, these dressed operators are used to calculate baryon matrix elements. Applications are worked 
out for the masses of the baryon octet, the meson-nucleon sigma terms, the magnetic moments of 
the baryon octet, the nucleon charge radii, the strong vector meson-nucleon couplings and the full 
momentum dependence of the electromagnetic form factors of the nucleon. 
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I. INTRODUCTION 



Chiral symmetry plays an important role in the low-energy (below 1 GeV) domain of Quantum Chromodynamics 
(QCD): it governs the strong interaction between hadrons. All known low-energy approaches (effective field theories, 
Lattice QCD, QCD sum rules, different types of quark models, etc.) in the study of the properties of light hadrons 
have to incorporate the concept of at least an approximate chiral symmetry to get reasonable agreement with data. 

The most convenient language for the treatment of light hadrons at small energies was elaborated in the context of 
Chiral Perturbation Theory (ChPT) [HQ, the effective low-energy theory of the strong interaction. ChPT is based 
on a chiral expansion of the QCD Green functions, i.e. an expansion in powers of the external hadron momenta 
and quark masses. It was proved that ChPT works perfectly in the meson sector (especially in the description 
of pion-pion interactions). A manifestly Lorcntz invariant form of baryon ChPT was suggested in Ref. In the 
baryon sector a new scale parameter associated with the nucleon mass shows up and this leads to certain difficulties 
in the formulation of a consistent chiral expansion of matrix elements. In particular, the chiral expansion of the loop 
diagrams starts at the same order as the tree-level graphs. This leads to an inconsistency in the perturbation theory: 
the higher order graphs contribute to the low-order ones and the physical quantities (e.g., nucleon mass) require 
renormalization at every order of the expansion. Later, a method, referred to as Heavy Baryon Chiral Perturbation 
Theory (HBChPT) 0, was suggested, which is able to avoid the problems with chiral power counting. HBChPT keeps 
track of power counting at every step of the calculation. A disadvantage of HBChPT is the lack of manifest Lorcntz 
covariance due to the nonrelativistic expansion of the nucleon propagator. As was pointed out in Ref. ^ this method 
also suffers from a further deficiency. Namely, the nonrelativistic expansion of the pion-nucleon scatteringamplitude 
generates a convergence problem of the perturbative series in part of the low-energy region. In Refs. [5|, |a, |3, l3] a new 
method for the study of baryons in ChPT was suggested. It is based on the infrared dimensional regularization of loop 
diagrams , which exploits the advantages of the two frameworks formulated in Ref. and Refs. 3 , while avoiding 
their disadvantages. An equivalent formulation of Baryon ChPT based on the extended on-mass-shell renormalization 
was suggested in Refs. 0, 0|. A successful application of the improved versions of Baryon ChPT to nucleon 
properties has been performed in Refs. [sl fol llOL llHIH. fTsf . In Ref. 0| the method has been extended to the multi- 
nucleon sector. The consistent inclusion of vector mesons in Baryon ChPT has been done in Ref. [3, which helped 
to successfully improve the momentum behavior of the nucleon form factors up to approximately 0.4 GeV^ [l5| . 

Unfortunately, in the context of Baryon ChPT one is able to calculate the momentum dependence of hadronic 
matrix elements only in a sufficiently narrow region (e. g., a n accurate description of nucleon electromagnetic form 
factors has been achieved up to about = 0.4 GeV^ |lOl lisjl'). Also, chiral symmetry is not the only important 
feature of strong interactions in the low-energy domain. There are the problems of hadronization and confinement 
which are completely avoided in the effective field theories dealing with hadronic degrees of freedom. These additional 
effects certainly have a strong impact on the hadronic interactions at intermediate energies. As an illustration of the 
importance of hadronization and confinement phenomena related to hadronic properties we recommend the review [T^ . 
In Ref. |l6j the general principles of QCD-motivated relativistic quark models have been formulated with the explicit 
incorporation of the three low-energy key phenomena: hadronization, confinement and approximate chiral symmetry. 
Many properties of light and heavy flavored baryons (including such sophisticated characteristics as slope parameters 
and form factors) have been successfully described in this approach and in a later developed similar model [r^ll8l |. 

The main objective of the present work is to develop a Lorentz covariant chiral quark model 19] which is consistent 
with the latest developments in the baryon sector of ChPT and leaves space for the additional features of low-energy 
QCD - hadronization and confinement. The full approach is an extension of the original chiral quark model 
suggested and developed in Refs. [20l l2ll |2^. We treat the constituent quarks as the intermediate degrees of 
freedom between the current quarks (building blocks of the QCD Lagrangian) and the hadrons (building blocks of 
ChPT). This concept dates back to the pioneering works of Refs. [23.l25j . Furthermore, our strategy in dressing the 
constituent quarks by a cloud of pseudoscalar mesons is motivated by the procedure pursued in Ref. [25| . Recent 
analyses of experiments at Jefferson Lab (TJLAB) 26], Fermilab BNL [23 and IHEP (Protvino) renewed the 
interest in the concept of constituent quarks. The obtained data can be interpreted in a picture, where the hadronic 
quasiparticle substructure is assumed to consist of constituent quarks with nontrivial form factors. These experiments 
also initiated new progress in the manifestation of constituent degrees of freedom in hadron phenomenology (see, e.g. 
Refs. Wl). 

The broader concept of chiral quark models dates back to the work of the early eighties 11^ 113, where the nucleon 
is described as a bound system of valence quarks with a surrounding pion cloud simulating the sea-quark contributions. 
These models aim to include the two main features of low-energy hadron structure, confinement and chiral symmetry. 
With respect to the treatment of the pion cloud these approaches fall essentially into two categories. The first type of 
chiral quark models assumes that the valence quark content dominates the nucleon, thereby treating pion contributions 
perturbatively [20I Isil Is^ . Originally, this idea was formulated in the context of the cloudy bag model I3H . By 
imposing chiral symmetry the MIT bag model |34| was extended to include the interaction of the confined quarks 
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with the pion fields on the bag surface. With the pion cloud treated as a perturbation on the basic features of the 
MIT b ag, p ionic effects generally improve the description of nucleon observables. Later, similar perturbative chiral 
models poTIs^ were developed where the rather unphysical sharp bag boundary is replaced by a finite surface thickness 
of the quark core. By introducing a static quark potential of general form, these quark models contain a set of free 
parameters characterizing the confinement (coupling strength) and/or the quark masses. The perturbative technique 
allows a fully quantized treatment of the pion field up to a given order in accuracy, usually evaluated in leading order. 
Although formulated on the quark level, where confinement is put in phenomenologically, perturbative chiral quark 
models are conceptually close to chiral perturbation theory on the hadron level. Alternatively, when the pion cloud 
is assumed to dominate the nucleon structure this effect has to be treated non-perturbatively. The non-perturbative 
approaches are based for example on Refs. |33|| . where the chiral quark soliton model was derived. This model is 
based on the concept that the QCD instanton vacuum is responsible for the spontaneous breaking of chiral symmetry, 
which in turn leads to an effective chiral Lagrangian at low energy as "derived" from QCD. On the phenomenological 
level the chiral quark soliton model tends to be advantageous in the description of the nucleon spin structure, that is 
for large momentum transfers, but is comparable to the original perturbative chiral quark models in the description 
of low-energy nucleon properties. 

As a further development of chiral quark models with a perturbative treatment of the pion cloud we 
extended the relativistic quark model suggested in for the study of the low-energy properties of the nucleon pH 
l23.l23l|. In the current manuscript we perform an extension of our previous approach in two directions: we formulate 
a chiral quark Lagrangian, which dynamically generates the dressing of the bare constituent quarks by meson degrees 
of freedom up to fourth order. We also formulate a manifestly Lorentz-covariant version concerning the structure 
of the bare constituent quarks, exemplified mainly for the case of the electromagnetic form factors of the nucleon. 
The resulting expectation values of dressed constituent quark operators evaluated for nucleon states are matched and 
checked in their low-energy behaviour with results of Baryon ChPT. 

In the manuscript we proceed as follows. First, in Section II, we derive an effective Lagrangian, which is taken from 
Baryon ChPT HHEIES, and formulate it in terms of quark and mesonic degrees of freedom by also including 
external fields. Second, we use this Lagrangian to perform a dressing of the constituent quarks by a cloud of light 
pseudoscalar mesons and by other heavy states. In this vein we use the calculational technique developed by Bcchcr 
and Leutwyler @ . We derive dressed transition operators with a proper chiral expansion, which in turn are relevant for 
the interaction of quarks with external fields in the presence of a virtual meson cloud. In Section III, we work out the 
implications for the chiral expansion of the nucleon mass and the meson-nucleon sigma terms, which, in consistency 
with Baryon ChPT, gives constraints on the parameters entering in the chiral quark Lagrangian. Applications to 
various quantities are worked out and presented in Section IV. We discuss the masses of the baryon octet and mainly 
the pion-nucleon sigma term. Using model-independent constraints on the bare constituent quark distributions in 
the octet baryons, we work out model predictions for the magnetic moments. Using finally a full parameterization 
of the bare constituent quark distributions in the nucleon, we give results for the full momentum dependence of the 
electromagnetic form factors of the nucleon and indicate the role of the meson cloud contributions. 

II. APPROACH 
A. Chiral Lagrangian 

The chiral quark Lagrangian Cqu (up to order p^), which dynamically generates the dressing of the constituent quarks 
by mesonic degrees of freedom, consists of the two main pieces Cq and Cjj'- 

Cqu = Cq + Cu, Cq = + c[^^ + cf^ + cf^ + . . . , Cu ^ C^^ + . . . . (1) 

The superscript (i) attached to ^q^u^f denotes the low energy dimension of the Lagrangian: 



= ^{u.u'^+X^), £(1) 
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^ ci {x+)qq - ^ {u^u,) {qD^^ D'' q + h.c. ) j- {u^u^)qq 



^ qi a>^- \u^, ^ q'^'^" 9 - qMq + Cr-,qx+q + . . . , (2) 



^ !^,-[^.^^+]^.,+ h.c. 
^ Zm 
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+ y (x+> ^^"^ 9 + 5 9 {^^+ x+} 9 + 1 9 'T^" X+) g - ^ g- [i?", [D^, F+]]a^'^ q+..., 

where x+ = X+ ~ i(x+) ? the symbols ( ), [ ] and { } occurring in Eq. |(2Jl denotes the trace over flavor matrices, 
commutator and anticommutator, respectively. 

The couplings m and g denote the quark mass and axial charge in the chiral limit, c,, di and are the second-, 
third- and fourth-order low-energy coupling constants, respectively, which encode the contributions of heavy states. 
Note, that the inclusion of higher-dimensional terms in the chiral quark Lagrangian was originally suggested in 
Ref. 25]. In particular, as a dimensional parameter in the higher-dimensional terms one can use the scale parameter 
of spontaneously broken chiral symmetry ~ AttF ~ 1 GeV {F is the octet decay constant 11, 36j) instead of the 
constituent quark mass. This replacement is equivalent to a redefinition of the values of the low-energy constants, 
e.g. C2 —> C2 (m/A^)^ , c?io — > dio{m/A^), etc. Both scale parameters m and A^^ are counted as the same order 
quantities in the chiral expansion, i.e as quantities of order 0(1). 

Additional terms in the Lagrangian denoted generically by dots do not contribute to the electromagnetic nucleon 
form factors and meson-nucleon sigma-terms, which are explicitly worked out later on in the applications. Here q is 
the quark field, the octet of pseudoscalar fields 

8 / n"/V2 + r]/VQ 7r+ K+ \ 

= ^0,A, = \/2 TT- -tt^/V^ + v/VQ • (3) 

r=i V K- K° -2r]/V6 J 

is contained in the SU(3 ) rn atrix U ~ ~ exp{i(l)/F) where F the octet decay constant 0,|23. We introduce the 
standard notations |3. Isl l35| 



lU' 



The fields i?^ and include sources 



Rf, = Vfj_ + a^, = eQ Afj_ - Q tanOwZ'^ + . . . , (5) 

= - a,, = eQA^ + ( . ^ A3 - q] tanOw ^" + — -^-=(W+T+ + h.c.) + . . . 

Vsm Ow / sm9w\2 

where s,p,v^ and denote the external scalar, pseudoscalar, vector and axial fields; 

Af^, and are the electromagnetic field, the weak charged and neutral boson fields. The quark charge matrix 
is denoted hy Q = diag{2/3, —1/3, —1/3} and 

Vud Vus 

T+ = I I (6) 


is the weak matrix containing the Cabibbo-Kobayashi-Maskawa quark- mixing matrix elements Vij . The tensor 
is defined as — u^Ffj^^Qu + uF^^Qu^ where F^i^ = d^A^ — d^Afj_ is the conventional photon field strength tensor. 
Here A4 = diagj™, tti, ms} is the mass matrix of current quarks (we restrict to the isospin symmetry limit with 
rhu — rhd ~ rh ~ 7 MeV and the mass of the strange quark rhs is related to the nonstrange one as rhs = 25 m). The 
quark vacuum condensate parameter is denoted by 

i? = -ij(0|S7.|G)--i^(0|Jd|0). (7) 

To distinguish between constituent and current quark masses we attach the symbol ' ("hat") when referring to the 
current quark masses. We rely on the standard picture of chiral symmetry breaking (B ^ F). In the leading order 
of the chiral expansion the masses of pseudoscalar mesons are given by 

= 2mB, Ml ^{m + rhs)B, = ^{m + 2m,)B. (8) 
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In the numerical analysis we will use: M^r = 139.57 MeV, Mk = 493.677 MeV (the charged pion and kaon masses), 
= 574.75 MeV and the canonical set of differentiated decay constants: = 92.4 MeV, Fk/F^^ = 1.22 and 
F^/F^ = 1.3 M 

The use of the physical masses and decay constants of pseudoscalar mesons incorporates only part of the corrections 
due to the breaking of unitary flavor symmetry. To generate another part of SU{3) symmetry-breaking corrections 
we added a string of terms to the Lagrangian Q): the current quark mass term qMq, terms containing LECs C5, 64, 
65, 67 and eg. The flavor-symmetry breaking terms containing term qA4q and LECS C5, 64 and 65 allow to decouple 
the mass of the strange quark from the isospin-averaged mass m. The fourth-order couplings ey and eg incorporate 
explicit SU{3) symmetry-breaking corrections in the magnetic moments of the constituent quarks and baryons. As 
was shown in Ref. [Tlj the inclusion of the SU (3) symmetry-breaking terms is sufficient to obtain agreement with the 
experimental data for the magnetic moments of the baryon octet. 



B. Inclusion of vector mesons 



Following Refs. 0, Hfl l38| we also include vector mesons in the chiral quark Lagrangian. In particular, we employ 
the tensor field representation of the spin-1 fields: vector mesons are written in terms of the antisymmetric tensor 
fields Wfii, = —W^^ where the three degrees of freedom Wij{i,j = 1,2,3) are frozen out. The latter representation 
is most convenient to construct the chirally invariant couplings of vector mesons to pions, photons and fermions 
(baryons, quarks): 

Cv = + C^' (9) 

where Cy is the free vector meson Lagrangian 

1 Af2 

£0, = --d^'W^^dpW'"'''' + ^W;^!^^"'" (10) 

and £y* = -f is the interaction Lagrangian of vector mesons with external vector and axial-vector sources 

(£y''^) and with baryons (£y*'^): 



^ -^pu j^^^ q + q^t^Spq + qa"" U^p [D^, [D'' , q]] . 

Here we define 



(pO+L^)/\/2 p+ K*+ 

W^.= \ p- {~p^+u)/V2 X*o I , (12) 

K*- K*° 



where Ri, Si and Ui are the effective couplings related to the ones {gvqq and ky) used in the canonical interaction 
Lagrangian of vector mesons with quarks (for details on the nucleon- level Lagrangian see Ref. [T(i | ) : 

Cy,,^'-^q{^rV,-^^a^^d.v)iq. (13) 

The standard nonet matrix of vector mesons is denoted by (its fiavor content coincides with the one of the 
antisymmetric tensor W^^)- The matching conditions relating the couplings are: 

i?S = 0, RT = -kygy,,^^, Us^-Ss, Ut^-(^^ + St). (14) 

4mV2 m m \My V2 / 

The couplings Fy and Gy are determined by the decays widths of p — > e'^e" and p —^ nn. Using low-energy theorems, 
e.g. p-nieson universality and the Kawarabayashi-Suzuki-Fayyazuddin-Riazuddin (KSFR) relation, one can express 
Fy and Gy by gygg and the vector meson mass My |33: Fy = My/gyqq and Gy — Fy/2. Hence, in the 
vector-meson sector we only deal with a single free parameter ky- 
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C. Power counting 

The Lagrangian set up for constituent quarks, pseudoscalar and vector mesons is assumed to be valid in the region 
between the chiral symmetry breaking A^^ and confinement Aqcd scales. Such a Lagrangian is non-renormalizable 
due to the existence of an infinite tower of higher dimensional terms. A solution to this problem can be achieved 
if there is a dimensional parameter in the theory which suppresses corresponding non-renormalizable terms in the 
matrix elements. The problem of power counting in non-renormalizable effective chiral quark theories was discussed 
in detail in Ref. [2^. It was shown in the framework of the so-called "naive dimensional analysis" that higher-order 
corrections in the matrix elements can be suppressed by a dimensional parameter A^ = AnF ~ 1 GeV. Note, that a 
deficiency of quark models is the smallness of the constituent quark mass. Formally, we need to treat this quantity as 
0{q'^), while the pseudoscalar meson masses as 0{q). Analytically the constituent quark mass survives in the chiral 
limit, while the pseudoscalar meson masses vanish. However, the numerical values of kaon and eta-meson masses 
are similar to the constituent quark mass. The contribution of the meson cloud is normally divided by the power 
of A;^ = AirFp where P = tt, K, rj. The constituent quark mass can be finally removed from the expressions for the 
observables using its universal relation to the nucleon mass which is valid in the chiral limit and at one-loop (see 
details in Sec. Ill): 

m rh \ g J 

In the last expression mjy and mjy are the nucleon masses at one-loop (the physical mass) and in the chiral limit, 
respectively, rh = to„ = is the nonstrange constituent quark mass at one-loop (the dressed nonstrange constituent 
quark mass). Finally, all expressions for the baryonic observables do not contain the scale parameter - the constituent 
quark mass, but contain the dimensionless parameter - the axial charge of the constituent quark. 



D. Dressing of quark operators 

The total effective Lagrangian C^g includes the two terms Cqu and Cy- The first term (Cqu) is responsible for 
the dressing of quarks by a cloud of pseudoscalar mesons and heavy states, while the second one (Cy) generates the 
coupling to vector mesons. Any bare quark operator (both one- and two-body) can be dressed in a straightforward 
manner by use of the effective chirally- invariant Lagrangian CcS- To illustrate the idea of such a dressing we consider 
the Fourier-transform of the electromagnetic quark operator: 



J,TnM - J d^x e-''- (x) , (x) = qix) 7, Q q{x) . (16) 

In Figs.l and 2 we display the tree and loop diagrams which contribute to the dressed electromagnetic operator J^'^om 
up to fourth order. 

The dressed quark operator jjj'cm(^) ^^'^ Fourier transform Jt^^cmi^) have the following forms 



q—u.d,s 



^^^d''[qix)a,M^)]} (17) 



<Tm('Z) = j d'xe-^'^-j^^X^ix)^ Jd^xe-^"^ E 



q—u.d,s 



7. /|,('Z') + ^^M. 9^/^(9') 



q{x) , 



where niq is the dressed constituent quark mass (see details in Sec. Ill); f]^{q^), ff){q^), fhil^) fp{<l^)i fpil^), 
fp{q^) are the Dirac and Pauli form factors of u, d and s quarks. Here we use the appropriate sub- and superscripts 
with a definite normalization of the set of /^(O) = (quark charges) due to charge conservation. Note, that the 
dressed quark operator satisfies the current conservation: 

d'^jtTlix)^ E [fU-d^)dn<i{^)lA^)] + ^-^^^^^ (18) 

Evaluation of the diagrams in Figs.l and 2 is based on the infrared dimensional regularization suggested in Ref. ^ to 
guarantee a straightforward connection between loop and chiral expansion in terms of quark masses and small external 
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momenta. We relegate the discussion of the calculational technique (T(l| to the Appendices A (infrared regularization) 
and B (exphcit form of the loop integrals). 

To calculate the electromagnetic form factors of the nucleon (or any baryon) we project the dressed quark operator 
between the nucleon (baryon) states. In the following we restrict to the case of the nucleon (the extension to any 
baryon is straightforward). The master formula is: 

W)\Jt':Z{q)\Nip)) = i2n)U\p' -p-q)rMp')\l,F,''{q') + a.^q^^ F,'' iq')]u^ip) 

I IrriN J 

= {2^f5\p'-p-q) J2 |/|.(9')(W)h-'r(0)l^(^'))+*|^/p('?')(W)|jX'(0)m^^^ (19) 

q—u,d 

where N(p) and ujsfip) are the nucleon state and spinor, respectively, normalized as 

{N{p')\Nip)) = 2EN{2nfS^p-p') (20) 

and 



u{p)u{p) = 2tojv (21) 

with Em being the nucleon energy Em = a/ to^ + p"^. Here {q^) and F^ (q"^) are the Dirac and Pauli nucleon form 
factors. For convenience we present the expressions for the isoscalar ^"^2) ^^'-^ isovector ^"^2) iiucleon form factors in 
Sec. IVc. In Eq. H19|) we express the matrix elements of the dressed quark operator by the matrix elements of the 
bare operators. In our case we deal with the bare quark operators of the vector jjj^g°(0) and tensor j^ffq{0) structures 
defined as 

J^r(O) = 9(0) 7m 9(0) , jX'(O) = 9(0) 9(0) • (22) 

Eq. (|19|l contains our main result: we perform a model- independent factorization of the effects of hadronization 
and confinement contained in the matrix elements of the bare quark operators ijj^q''(0) and jJj^'^J(O) and the effects 
dictated by chiral symmetry (or chiral dynamics) which are encoded in the relativistic form factors f%{q^) and fp{q^). 
Due to this factorization the calculation of f1){q'^) and fp{q^), on one side, and the matrix elements of jjj^^''°(0) and 
j^ffq{0), on the other side, can be done independently. In particular, in a first step we derived a model-independent 
formalism based on the ChPT Lagrangian, which is formulated in terms of constituent quark degrees of freedom, for 
the calculation of /|)((?^) and fp{q^)- The calculation of the matrix elements of the bare quark operators can then be 
relegated to quark models based on specific assumptions about hadronization and confinement. The explicit forms of 
f1){q^) and fp{q^) are given in Appendix C. Note that we also show in Fig. 3 the diagrams contributing to the strong 
vector meson- nucleon interactions pNN and uNN at one loop which will be discussed in Sec. IV(D). 



E. Matching to ChPT 



The matrix elements of the bare quark operators should be calculated using specific model assumptions about 
hadronization and confinement. However, the use of certain symmetry constraints, discussed in the following, leads 
to a set of relationships between the nucleon and the corresponding quark form factors at their normalization point 
at zero momentum. In general, due to Lorentz and gauge invariance, the matrix elements in Eq. H19|) can be written 
as 

{Nip')\jlT{0)\N{p)) = SAr(p')(7M^f'(9') + :^^,.'l''F,^Hq')]uN{p), (23) 
^-^(W)|JX°(0)l^(p)) = ^Ar(p')(7MGf'(9') + :r^a,,q-Gp{q')]uMip), 

where f/J^' (g^) and Gf[^)(g2) are the Pauli and Dirac form factors describing the distribution of quarks of flavor 
q = u,d in the nucleon. 

The first set of relations arise from charge conservation and isospin invariance: 

Ff"(0) = Fl"^{0) = 2 , ^['^(O) = Fr(0) = 1 , Gf «(0) = , (24) 
F|"(0) = i^2"'(0), Ftm^FnO), Gr(0) = Gf(0), Gf{0)^GT{0). 



8 



Note, that the quantities G^'^(O) are related to the bare nucleon tensor charges S^'^'^'^- 



where (5^^™ are defined by ^0 



GT{Q) = (0) = ^5^r = "^SlT. (25) 

Gf (0) = Gr(o) = — <5^r = — C' 

m ' m 



The second set of constraints are the so-called chiral symmetry constraints. They are dictated by the infrared-singular 
structure of QCD in order to reproduce the leading nonanalytic (LNA) contributions to the magnetic moments and 
the charge and magnetic radii of nucleons [lfiH4ll|: 

/„2nM ^ 9a 

where gA and tojv are the axial charge and mass of the nucleon in the chiral limit. In particular, the LNA contribution 
to the magnetic moments is proportional to M^^. The nucleon radii are divergent in the chiral limit. The LNA 
contribution to the charge radii is proportional to the chiral logarithm ln(M7r/ mpf). The LNA contributions to 
the magnetic radii are represented by the same logarithm as in the case of the charge radii and by the singular 
term proportional to 1/M„. In order to fulfill the chiral symmetry constraints (|27|) we derive the following identities 
involving the F^'(O) and G^'(O) form factors and the low-energy constant (LEG) dio: 

1 + FP«(o) - Ft{0) = Gf (0) - Gf{0) ^(^] — , (28) 

\ 9 J m 

1 + F^'^iO) - F2""(0) = Gf{0) - Gr(0) 



2 



and 



9 J m 

Here d^^^'^ and Jio denote the renormalized LECs dg^^""" and dio, respectively, at /i =mN: 

^-ChPT^^ChPT^l+M^^ ^ + ii^A (30) 



^(m) = -^(t^ - i(ln47r + r'(l) + 1)L A = A(m^.) . (31) 



where 

(47r)2 {d-A 2 

Identity H29|) represents the matching condition between the LEG of the ChPT Lagrangian and our quark-level 
Lagrangian. Analogous conditions involving other LEGs can be derived when matching other physical ampli- 
tudes/quantities (see, e.g. Sec. IIII|I . 

Applying SU(6)-symmetry relations of the naive nonrelativistic quark model for the ratios of the magnetic moments 
and the tensor charges of the nucleon we can derive additional and well-known relations between F^*(0) and G^'(O), 
respectively: 

2 + Fr(0) _ 2-fF,"'^(0) _ Gr(0) _ Gno) _ ^ , . 

i + Ff(o) i + Fr(o) Gf(o) Gr(o) ■ ^ ' 
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Substituting Eq. ll^ into Eq. if^ . we arrive at 



4 / V 



5V .9 
1 / \ ^ mw 



5V 5 



5 V 5 / 

Gr(0) = Gf(0) _ 

GfiO) = Gr(0) = - i (^y ^ ' ^ . (33) 
Other interesting results are the predictions for the bare values of nucleon magnetic moments and tensor charges: 

Mr'^ ^ -l^^^n"' = E e.[^r(o) + F^m = U^-fX^ (34) 

q=u,d \ 9 / 

and 

<s^^r ^ -45pY'' = = t (7)'' (35) 

where e„ = 2/3 and — —1/3 are the electric charges of u and d quarks. Below, in Sec. IIIII we also derive the 
constraint relating the axial charge of the nucleon 5^ and of the constituent quark g. 

III. NUCLEON MASS AND MESON-NUCLEON a-TERMS 

In this chapter we consider two other important quantities of low-energy nucleon physics: the nucleon mass and 
meson- nucleon cr-terms which are constrained by the Feynman-Hellmann theorem (FHT) 42]. In particular, the FHT 
relates the derivative of the nucleon mass with respect to the current quark masses to the pion-nucleon sigma-term 
cr^AT and to the strange quark condensate in the nucleon: 

<J^N un{p)un{p) = rh{N{p)\u{Q)u{Q) + d{Q)d{Q)\N {p)) = um{p) un{p) , (36) 

om 

y,UN{p)uN{p) = {N{p)\s{0)smN{p)) = ^^UN{p)uN{p). 

In quantum field theory the nucleon mass mjv is defined as the matrix element of the trace of the energy-momentum 
tensor Q^'''{x)■. 

mNUN{p)uN{p) = {N{p)\Q^;{0)\N{p)) . (37) 

When we restrict to one-body interactions between the quarks and neglect to the contribution of the heavy quarks in 
the constituent quark (CQ) approach the master formula H37(l spells as 

mjv un{p)un{p) = (iV(p)|7^mas«(0)|A^(p)) . (38) 

'^mass(a;) = q{x)mqq{x) is the part of the Hamiltonian referred to as the quark mass term where niq — 
diag{m,i, TO(j, m^} is the matrix of constituent quark masses with m„ — md = fh due to isospin invariance. Note, 
Eq. (|38|l is valid in the one-body approximation. In general, one should include in the trace of the energy- momentum 
tensor also two- and three-body quark operators which encode quark-quark interaction contributions to the nucleon 
mass. In some CQ models the nontrivial dependence of ruq (or rriAr) on the current quark masses is missing. This 
leads to a contradiction with the low-energy behavior of the nucleon mass as a function of mirhs) and with the 
Feynman-Hellmann theorem. The use of the effective Lagrangian constrained by Baryon ChPT enables one to 
perform an accurate and consistent calculation of the nucleon mass and the corresponding sigma-terms. Again, as in 
the case of the electromagnetic form factors, the extension to other baryons is straightforward. 



10 



In analogy with the electromagnetic form factors we define the nucleon mass and later on the sigma-terms as 
expectation values of the dressed operators. First, we write down the bare quark mass term: 

nt^l!.ix)^mq{x)qix) (39) 

i.e. the quark mass term at leading order of the chiral expansion (in the chiral limit). Here m is the value of the 
constituent quark mass in the chiral limit introduced before in the Lagrangian The nucleon mass in the chiral 

limit mjv is defined by 

mN um{p)um{p) = {N{p)\-H'::^ltM\N{p)) = m{N {p)\m {p)) ■ (40) 

The dressed quark mass term and the physical nucleon mass are given by 

<all(2;) = qix)m,q{x), (41) 
mNUNip)uN{p) = (7V(p)K,^||(0)|iV(p)) ={Nip)\mm,qmN{p)) 

where uiq = diag{m„, m^, TOs} is the matrix of the dressed (physical) constituent quark masses (in our case the 
constituent quark mass at one loop with inclusion of chiral corrections) with m„ — rrid due to isospin invariance. 
The constituent quark masses ruq = mq{rh,rhs) have a nontrivial dependence on the current quark masses m and rhs 
which can be accurately calculated with the use of the chiral Lagrangian For illustration we discuss the explicit 
expressions for the nonstrange and strange constituent quark masses at one loop and at O(p^) with 

rUq = m + Tiqlm) . (42) 

The quark mass operator = diagjEu, S^, Ss}, with = E^ = 5] due to isospin invariance, is evaluated on 
the mass-shell p = m, which is ultraviolet-finite by construction. All ultraviolet divergencies are removed via the 
renormalization of the fourth-order LECs ei, 63, 64 and 65 contributing to the E^ operator. Here and in the following 
we identify the quark mass occurring in the loop integrals with its leading order value niq — > m. The operator Eg is 
described by the diagrams in Fig. 4 and after expansion in powers of meson masses is given by: 



(43) 



and 



327r \ 3 F|. 9F2 J 647r2m\ 3 F|. 9F2 

o 3c2 f Ml 4 Mi M^ 14,, 

^^^^ + 1^1^ + 3^ + 3^h r^^''- - ''^^ 

e-.M^ + |(M^ - Mlf - ^-lM\Ml - M^ + |(A4 - Af^)^ 



E _ ^ _ V r4M| + 4M|l ^r4< 4<1 

^ " ^ 32^ \3 Fl 9 F^ j 64^2^13 ^2 + 9 ^2 / ^^^> 

Ac Af2 , r A/4 4 Ml M^ 1 8 2 ^,2^ 

- ^'^^^ +128^17^ + 3^1 +37|rr^^^^^~^-^^) 

+ e,M^ + |(A/i - A/J)2 + ^MHmI - Ml) + |(A4 - Mlf 
6 3 9 



where 



M^ = Ml - Ml + ^M^ ^ ^Ml + Ml . (45) 

Note, E and Eg degenerate for m = rhs and Af2 coincides with Ml at TOs = 0. The contribution of the diagram in 
Fig. 4(5) (the so-called mass-insertion term) is generated by the terms from second-order Lagrangian and is absent 
in final results (see Eas. (|43|) and 144|l ') after certain replacement the constituent quark mass in the free Lagrangian 
(see details in Q). For simplicity the chiral logarithms are hidden in the renormalized LECs e,; with i = 1, 3, 4, 5: 

= - 32^^"^ ('^^ 
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where 

A.^(-^{^-^(ln4. + r'(l) + l)|. (47) 

The LECs contain the poles at d = 4 which are cancelled by the divergencies proportional to in the r.h.s. of 
Eq. H46() . Therefore, the renormalized couplings e[(/i) (or e^) are finite. The set of the coefficients is given by 

R 205" 32 4g2 16 

= yy - y^^ ' = -yr - y^^ ■ (^^^ 

Note, that the — m splitting is mainly generated by the difference in the values of the strange and the nonstrange 
current quark masses. The chiral symmetry constraints and the matching of the nucleon mass calculated within 
our approach to the model-independent derivation of Q allow to deduce certain relations between the set of our 
parameters and the ones in Baryon ChPT. In particular, the coefficient connected with M^/F^, referred to in the 
literature as the leading nonanalytic coefficient (LNAC) is model-independent and constant due to dimensional 

arguments. Precisely, this coefficient is equal to —3g\/{327rF^). The cubic term in powers of the meson mass shows 
up due to the infrared singularity of the diagram in Fig. 4(3). Using this requirement, we can relate the value of axial 
charge of the constituent quark in the chiral limit to the corresponding nucleon quantity: 

gluN{p)uN{p) = 9^N{pM0)qmN{p)) ■ (49) 

The matching condition gives a constraint for the matrix element of the bare scalar-density operator in the nucleon. 
A rough estimate with gA = 1.25 .3|, 5 ~ 1 and by taking into account the normalization of the nucleon spinors gives 
a quite reasonable value for the scalar condensate: {N{p)\q{0) q{0)\N{p)) ~ (25/8) toat . Using the matching condition 
()49|l we derive a final expression for the nucleon mass at one loop: 



where 



and 



O 

TTlAT —TUN (50) 



mjy — [ — m, mj^ ~ — \ m (51) 

9 J \ 9 J 



— oo_ 1 r2 + o P2 + n 172 f G/i_2™ 1 i?2 + o i?2 + n i?2 I V^^j 





35i 


2M4 




9F^j 


647r2TO \ F2 


3Fi 


+ 


3C2 






4 


1287r2 ^ 


[F^ 3 Fl + 




3 



327r i F2 3 F2. 9F2 j 647r2TO { F^ 3 F^ 9 F^ 

s 2 r 

9A 
9 



m - 4ciM2 + — ^ ^ + ^ + ^ + c5(Af^ - Af2) 



+ e-iAf^ + |(A4 - Mlf - '-^M^Ml - Ml) + |(Afi - M^f 

The constituent quark mass can be removed from the expressions for nucleon observables using the matching condi- 
tions (|51|) . The same is true also for other baryonic observables and the constituent quark can be removed from the 
corresponding expressions. In the SU(2) picture (when neglecting the kaon and 77-meson loops and putting C5 and 65 
to be equal to zero) we reproduce the result of ChPT for the nucleon mass at one loop and at order 0{p*) Q 



mN = 



-Ac?^^^ - If H + hM^^ + k,Mt + OiM^) , (53) 



I. _ 3 / Of ChPT , ChPT I ^„ChPT\ 

- 32.2fH™. ^ ) 

I. _ ^ChPT 3 / 2g\ ChPT\ 



-ChPT _ gChPT _ 3A / g\ ^ SfChPT _^_ ^ChPT _^_ 4f ChPT^ 



4c^ 
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if we fulfill the following matching conditions between the LECs of the ChPT Lagrangian and our quark-level La- 
grangian: 



2 

9 A ' 



o ChPT ChPT /f^ChPT 9a _ 

aci — C2 — - — = 

-ChPT 3 / 2g\ ^ ChPT\ _ 



Ci - C2 - 4C3 

3 



a: 


( 


mN- 


\ 9 






o 









(54) 



2 

9A 



' 64^2^2 

where e* is a SU{2) analogue of ei derived in Eqs. (|46|l and H48|) in the three-flavor picture: 

el^el-^X^, = ^(^^-8ci+C2+4c3^ . (55) 

Note, that there is an additional condition on g and gA which shows up in the calculation of the axial nuclcon charge: 

gA um{p) 7^ 7' ^3 um[p) - g{N{p)\m 7^ 7^ T3 <7(0)|7V(p)) . (56) 

Eq. H56() gives a constraint on the matrix clement of the isovector axial current. 

As an example for the sigma-terms we consider the pion-nucleon sigma-term cr^jv- In QCD (see Eq. (|36|l ') this 
quantity is related to the expectation value of the scalar density operator. It is connected to the derivative of the 
part of the QCD Hamiltonian, explicitly breaking chiral symmetry, with respect to the current quark mass. This 
definition is consistent with the FH theorem ji^. In the context of CQ models, we should proceed with the dressed 
Hamiltonian 'H'^^Hix) which already showed up in the calculation of the physical nucleon mass. In particular, the 
dressed scalar density operators if^°^^{x) (where i — u,d,s is the flavor index), relevant for the calculation of the 
meson-baryon sigma-terms within our approach, are defined as the partial derivatives of Timass (^) with respect to the 
current quark mass rhi of z-th flavor: 

_^.dre.s(^) = ^Z1|PM ^ ,-(^) ^ ,(^.) ^ ,-(^) ^ q(^) . (57) 

orrii arm orrii 

In the case of ajrN we have: 

a.N un{p) UNip) = m (7V(p)| ^'^--(0) -t- <=-(0) |7V(p)) . (58) 
It should be clear that Eq. (|58|l is consistent with the FH theorem: 

T.N UMip) um{p) = ^ ^ ^(iV(p)K-(0)|jV(p)) = ^ ^ UNip) un{p) . (59) 

= mjv Sjv(p) tijv(p) 

Below we give the definitions of the strangeness content of the nucleon y^, the kaon- nucleon axN and the eta- nucleon 
ariN sigma-terms: 

^ dmN/drhs 
= 2 



drriM /dm 
rh + rris 



ur,{p)u^{p)a-^''^ = -1^{N{p)\j1\7)%+{^)\N{p)) , 
UNip)uN{p)air^ = !^^(iV(p)btr?(0) + jtr+(0)|iV(p)) 





2 


m 


+ rris 




4 


m 


+ rhs 


2 



UNip)u^{p)ajri^ = —^{Nip)\CJ:^^^iO)\Nip)), 

un{p)un{p)<J,n = i(iV(p)|mj^;-«(0) -f 4m,j^-(0)|iV(p)), (60) 

where = j^™**^ ± and \N{p)) denotes a proton state. All further details can be found in our previous 

papers |22| and in Appendix D. 
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IV. PHYSICAL APPLICATIONS 



In this chapter we consider the appUcation of our approach to specific nucleon properties: magnetic moments, 
electromagnetic radii and form factors, nucleon mass, meson-nucleon sigma-terms and strong vector meson-nucleon 
form factors. We also calculate some canonical properties of hyperons: masses and magnetic moments. Note, that 
these properties have been studied within different approaches: QCD sum rules kj, ChPT and HBChPT ^,^,^^01 
m El m mill 113, lattice QCD "isl, approaches based on vector dominance 1491 and on solutions of Schwinger- 
Dyson and Bethe-Salpeter equations |^] , different types of chiral and quark models 0, 0, |^ |S 1^1 El IS 1^ , 
etc. Our main interest will be restricted to the region of small external momenta where the contribution of the meson 
cloud is supposed to be largest. 



A. Masses of the baryon octet and the meson-nucleon sigma-terms 

The nucleon mass and the meson-nucleon sigma-terms at one loop depend on the set of parameters: g, m, ci, C2, C5, 
ei, 63, 64 and 65. First, we discuss a choice for g and m. Note, that these parameters are constrained in our approach 
by the matching condition H51|l : rriN — (gA/ gY' fn- In literature the value of the axial charge of the constituent quark 
varies approximately from 0.9 to 1 (see detailed discussion in Refs. |52j). The nucleon mass in the chiral limit was 

estimated in HBChPT: niN^ 770 ± 110 MeV ]^ and mN= 890 ± 180 MeV li^. Inserting this range of values for g 

and niN and also the experimental value for gA = 1.267 l53'| into Eq. H51(l gives the following limits for the constituent 
quark mass in the chiral limit: 

m ~ 500 ± 167 MeV . (61) 

Finally, we choose: 

TO = 420 MeV and g = 0.9. (62) 

Using the experimental value for the nucleon axial charge gA — 1.267 [sl and Eq. we get mN= 832.4 MeV 
which is rather close to our previous estimate done in the framework of the perturbative chiral quark model (PCQM): 
TOjv= 828.5 MeV Here we do not pretend on a more accurate determination of m and g. This will be done 

in a possible forthcoming paper where we intend to involve more constraints from data. However, we stress that we 
need a rather small value for g in the interval 0.9 — 1 to justify perturbation theory. In particular, the use of g = 0.9 
gives the shift of the nonstrange constituent quark mass equal to Am = 53 MeV. For g — 0.95 and g = 1 we get 
Am = 107 MeV and Am = 164 MeV, respectively. On the other hand, the choice of g is constrained by the bare 
magnetic moments of the nucleon (see Eq. ()76f) '). The use oi g = 0.9 gives a reasonable contribution of the valence 
quark to the nucleon magnetic moments (see discussion in the next section). Finally, the couplings ci, C2, C5 and ei 
are fixed by four conditions: 

TOAT = 938.27 MeV, cr^Ar=45MeV, = 0.2 , to^ -to -170 MeV. (63) 

The result of the fit is: 

Cl = (-0.317- 0.004^3 + 0.020^4 + 0.001^5) GeY~\ 
C2 = (1.093 - 1.354 §3 -I- 0.474 §4 - 0.338^5) GeV^S 
C5 = (-0.316 + 0.063^3 + 0.005^5) GeY^\ 

ei = (1.157 + 0.298 §3 + 0.162 §4 + 0.088 65) GeV~^ (64) 

where for convenience we introduce the dimensionless parameters = x 1 GeV"^. Note, the mass difference 
TOs — TO„((i) is a crucial quantity to roughly describe the splittings between the octet baryon masses. Latter quantities 
are proportional to TOs — to. Neglecting isospin-breaking and hyperfine-splitting effects we have in our approach: 

TOA - toat = - to ~ 170 MeV (data : 177 MeV) , 

TOE - TOAT = TO^ - TO ~ 170 MeV (data : 251 MeV) , (65) 

TOH - TOAT = 2[tos - to] ~ 340 MeV (data : 383 MeV) . 

The meson cloud corrections are important to reproduce the full empirical value of cTt^n, contributing about 2/3 to 
the total value (for a detailed discussion see Ref. ^3|)- In particular, our result for aj^M is compiled as: the total 
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value is a-^N = 45 MeV, the contribution of the valence quarks is a™^ = 13.87 MeV (30% of the total value), the 
contribution of the pion cloud is dominant 

<jv = 26.69 MeV + 0.45 MeV ^3 - 2.58 MeV 64 - 0.11 MeV 65 , (66) 

the kaon and 77-meson cloud generates 

cr^+'J = 4.44 MeV - 0.45 MeV + 2.58 MeV ^4 + 0.11 MeV 65 . (67) 

The separate contributions a'^j^ and t^j^^ depend on the parameters 63, 64 and 65, while the total contribution is 
independent on with i = 3,4,5. Our predictions for the kaon-nucleon axN and eta-nucleon ar/N sigma-terms are: 
cr]^^ = 381.9 MeV, cr^^ = 351.8 MeV, ajf^ = 366.8 MeV, crj^^ = 15 MeV, and a^N = 90 MeV. 



B. Magnetic moments of baryon octet 



The magnetic moments of the nucleon can be written in terms of the Dirac and Pauli form factors as 

MJV = <(0) + F2~(0), (68) 
where in our formalism F^{0) and i^2^(0) are of the form 

Ff{0) - E /|,(0)Ff''(0), 



q—u,d 



i^2~(0) = ^[/|,(0)Ff^(0) + /^(0)G^^(0)] 

q—u,d 



(69) 



where /|)(0) = Cq is the quark charge due to the charge conservation. It is convenient to separate the expressions for 
(0) and {0) into two contributions: the bare (valence-quark) and the meson cloud (sea-quark) contributions. 
The bare contribution is 



bare 



t^A^ bare 
^2 



(0) 
(0) 



q—u.d 



(70) 



The meson cloud gives rise to 



F, 



TV cloud /n^ — 



71 A'^ cloud 



(0) ^ 0, 

(0) ^ E /p(0)G2^H0) 

q—u^d 



(71) 



Note, that the contribution of the meson cloud to the Dirac form factor ^4^(0) is zero due to the charge conservation 
of the nucleon. Therefore, the nucleon magnetic moments are given in additive form with 



where 



,, bare , cloud 

MAT - Mat + Mat • 



,,barc J7A^bare/n^ i ^A^bare/n^ 

Mat = ^1 W+^l (U) 



(72) 
(73) 



and 



cloud 
Mat 



N cloud 



(0) + K 



Af cloud 



(0) 



(74) 



With the constraints laid out in Sec. II, we derived a relation for the bare contributions to the nucleon magnetic 
moments [see Eq. (I34|l ]: 



bare _ _ bare 
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(75) 
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Note that Eq. (|75|l can be further simpHfied using the constraints of (|51|l : 

^W^_3^ba.e^ 3j^^y^ (76) 

As aheady mentioned in the previous section, the choice of g is constrained by the bare magnetic moments of the 
nucleon (see Eq. H76|l ). The value of g = 0.9 results in a reasonable contribution of the valence quarks to the nucleon 
magnetic moments with: 



Mp^''^ M^" = ^ (^^ V ^ 2.357 (77) 



which is about 84% (for the proton) and about 82% (for the neutron) of the experimental (total) values, where the 
remainder comes from the meson cloud. 

The main parameters contained in the meson cloud piece that play an important role in fitting the magnetic 
moments are the second-order coupling cg and the fourth-order flavor- breaking LECs 67, and eg. The coupling eg is 
absorbed in cg after the certain redefinition. To constrain these values we need besides the experimental values for 
the nucleon (p and n) an additional value from the baryon octet. Therefore, we will proceed with an extension of our 
formalism to calculate the magnetic moments of the whole baryon octet. 

In the following, for convenience, we introduce the bare Dirac {Fi , Gf ) and Pauli (Fj , G^) form factors of the quark 
of flavor q: 

(^bOlj'rWkb)) = u,ip')i^^,F?iq^) + -l-a,,q'^F^iq^)y,{p), (78) 
*l^(g(y)|jX°(0)kb)) = u,ip')(^,Gl{q') + -^a,.q''GUq')]u,ip). 

Note that the valence quark form factors (or valence quark contributions) are constrained by certain symmetries 
including the infrared singularities (see discussion in Section II. E). Introducing the valence quark form factors we 
fulfill matching conditions between our approach and ChPT and, therefore, the structures due to chiral symmetry are 
not violated. 

The Sachs form factors of the quark of flavor q are: 

•^f (0 = F^'^it) + Gf (t) , .Ff (t) = Ff (t) + Gf (t) , (79) 

where 

F{G}f (i) = F{Gm - J-^F{Gm , F{G}f (t) = F{G}?(t) + F{G}l(t) , t = -q\ (80) 

are the contributions to the Sachs form factors associated with the expectation values of the vector and tensor currents, 
respectively. 

By using SU(6)-symmetry relations one can relate the Dirac and Pauli form factors describing the distribution of 
quarks of flavor q ^ u,d,s in the baryon "S", that is F^^^{t) and G^2)(i)j to Fq''^^\t) and Gq'"'^^\t) by 



pBq 



'F^{r)+a1?x'"'FtHi)rBY 



it) = Y^{-«|''Ff(t)+af/x^^Ff(t)^ 
Gf^(i) = ^^|a|«Gf(t)-faf/x^''Gf(i)rB 
G^'it) = -J_|-«|''Gf(t)+af/x^^Gf(t)|, (81) 



where tb = t/(4m^) and tub is the baryon mass. In addition to the strict evaluation of SU(6) we have introduced 
the additional parameter x^'^ for each quark of flavor q. The interpretation for adding these factors is such that to 
allow the quark distributions for hyperons to be different from that for the nucleons. In the case of the nucleons we 
set x^'^ = 1- The values for a^'^ and a^/ for the baryon octet as derived from SU(6)-symmetry relations are given in 
Table 1. 
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The quark Sachs form factors are modeled by the dipole characteristics with damping functions of an exponential 
form. This phenomenological form is required to reproduce in particular the deviation of the electromagnetic form 
factors of the nucleon from the dipole fit as evident from recent experimental measurements. We use the parameteri- 
zation 

- ^^^^^ - Pg'(^) 

0?W^7,<«)p^^. GfW = .?|^^||-p, (82) 

where Pq{t) = exp{—t/X'^^) and Pq'{t) = exp(— t/A^^j). The parameters /i^ and /i^ are fixed by the symmetry 
constraints [see Eqs. and lf75)) ]: 

A^,"^ - - Mp"°- (83) 

The remaining parameters 7^, Aq£(jv,/) and XqE{M) ^re to be fixed later when we consider the full momentum depen- 
dence of the nucleon electromagnetic form factors. Note, that in Ref. 54] a similar parametrization of the nucleon 
form factors has been considered. In Ref. [ssf the damping functions p{t) have been parametrized with constant 
values. 

The magnetic moment of the octet baryon " _B" can be written in complete analogy to the nucleon case as 

,, bare I cloud tiiA\ 

where 

q—u.d,s 

^cioud ^ J2 fUO)G^\0). (85) 

q—u.d,s 

At i = 0, the parameters cg, eg, 67, eg, C2, and ky contained in /p(0) will contribute, where the renormalized LECs 
eg, 67, eg are given in Appendix C. The parameter eg can be absorbed in the coupling cg via 

cg ^ cg = cg — IQmM^eQ . (86) 

Finally, form factors fp{0) depend on the parameters cg, 67, eg, C2 and ky- Three of them (cg, 67, eg) can now be 
fixed using the experimental values for the magnetic moments of the nucleons and of a hyperon in the baryon octet. 
We choose the A°-hyperon with /iao = —0.613 ± 0.004 (in units of the nuclear magneton). In our analysis of the 
magnetic moments of the octet baryons we present two cases. First, we restrict to the SU(6) case by setting all values 
X^* — 1. Secondly, we allow the to be additional parameters. From now on we will name these two cases shortly 
as "Set I" and "Set 11", respectively. 

A variation in the value of the axial charge g of the constituent quark will also give rise to different contributions 
of the bare and the meson cloud parts to the total values of the magnetic moments. However, the total magnetic 
moments are the same in any case. As was discussed in the previous section, for the axial charge we choose g — 0.9. 
In case of Set I we then obtain for cg, 67, and eg: 

£g = 0.859 - 0.106 52 - 1.254 fcv, 

67 = (-0.656 -I- 0.073 C2 -I- 0.607 fcy) GeV"^ (87) 
eg = (0.033 - 0.008 52 - 0.141 fcy) GeV-^ 

where 52 = 02(1 -I- fcy). The resulting values for the magnetic moments of the baryon octet for this case (Set I) are 
shown in Table 2, where reasonable agreement with data is obtained. Meson cloud contributions to the total values 
of the magnetic moments are about 5 — 30% depending on the baryon. 

For the case of Set II we start from the SU(6) result with the additional breaking parameter for fJ^'^'^ which 
can be written in terms of Pp^'^° as /ij^o"^ — — /ip^''°x^'*/3- Then we have 

,,barc 

X^^=-3^«-3^, 
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where we further assume that the contributions from the meson cloud part to the total magnetic moments both of 
A*^ and of p are of the same order. With the experimental values for /i^o and fip we get = 0.66. By fitting /ip, 
fin, and /xao, with = 0.66 as an additional input, we obtain 

cq = 0.834-0.106 52 - 1-254 fcy, 

67 = (-0.832 + 0.073 £2 + 0.607 fcy) GeV-^ (89) 
eg = (0.051 - 0.008 22- 0.141 fcy) GeV-^ 

where the constants contained in these parameters are slightly different from Set I. With appropriate values for the 
remaining , the central values of the experimental results of the other magnetic moments in the baryon octet can 
be fully reproduced. For the set of x^'^ we get 

X"" = x""* = 0.633, X"' = 0.694, 
^SA« ^ x^'^'' = 0.988. (90) 

relying on isospin symmetry. The isotriplet S+, and shares the same set of x^' for the quark of flavor g, while 
'EP and contains the same set of x"^- The parameters x^^" and x^^'^ are directly related to the E — A magnetic 
transition moment. For completeness, our corresponding results for the magnetic moments are also listed in Table 2. 
In the following we will use the last set of values for ce, ey, and eg, when considering the electromagnetic nucleon 
form factors. 



C. Nucleon electromagnetic form factors 

Using the parameters from the analysis of the magnetic moments and the masses of the baryon octet, further 
supplied by the ones from the meson-nucleon sigma-terms, we now can consider the full nucleon electromagnetic form 
factors. We recall that the Dirac and Pauli form factors for the nucleons are 

q—u^d 

F^{.t)= Y.^^D{i)F^''it) + mG2'{t)\. (91) 

q—u,d 

where F^^^{t) and G^2)(0 refer to the bare constituent quark structure, while /^(<) and fp{t) contain the chiral 
dynamics due to dressing of the quark operators. The isoscalar Ff{t) and isovector Ff{t) nucleon form factors are 
written as 

Ffit) = Ff(i)+F«(t)= E [fUt)iFrit)+Frit))+fUt)iGrit)+Grm (92) 

q—u,d 

FYii) = Ff{t)~Fni)= E [/^w(^rw-^rw)+/^w(Grw-Grw)] 

q—u,d 

where z 1, 2. For clarity we recall the forms of F^^^{t) and Gf^^^{t) of Eq. (jSlfl with: 

F^'it) = [^%'F^ + <Ft^r^'^ , F^^it) = ^ l^-a^^^Ff + } , 

Gi\t) = YTT;^ {^e'G^ + «A/Gf r^.} , G^\t) = ^ + a^'G^^^ , (93) 

where a^' and a^'' are the SU(6) spin-flavor couplings {a^ — — 2, a^^ — — 1, a^" — a^f — 4/3, 
a^^ = a^f = —1/3); t — —q^ is the Euclidean momentum squared and tat = t/(Am?^). The Dirac and Pauli form 
factors f1)(^p-jit) of the quark of flavor i are shown explicitly in Appendix C. Again, Eq. (|93|l can be separated into a 
bare part and a meson cloud part as 

q—u,d 

^Arba.c(^) = E ^^'^F,^'>it), (94) 

q—u.d 
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and the meson cloud contribution is 

= ^ [(/|,(i)-e,)f^f'W + /^(i)Gf'(t)], 

q—u.d 
q—u,d 

where are the electric quark charges. 

In the following we try to achieve a reasonable description of the electromagnetic form factors of the nuclcon by 
an appropriate set of parameters contained in F(^{t) and F^ (t). In general both the bare and the meson cloud parts 
will contribute to the form factors. The bare part should well describe the form factors up to high whereas the 
meson cloud part should play a role only for small Q^. Our strategy is that we first fit the bare part of all the form 
factors to the experimental data from ~ 0.5 GeV^ to high Q^. After that we reconsider the low region by 
taking into account the meson cloud effect. At this point we do not pretend on an accurate analysis (or prediction) 
of the nucleon form factors at large Q^. This is a task more appropriate for perturbative QC D (pQC D) which is the 
most convenient theoretical tool in this momentum region (for recent progress see e.g. Refs. pA \^^). In particular, 
in Ref. [H^l the following large Q^-behavior for the ratio of the nucleon Pauli F2 and Dirac Fi form factors has been 
derived using pQCD: 

F, In^(QVA^) 

where A is a soft scale related to the size of the nucleon. Our main idea is to use a physically reasonable parametrization 
of the nucleon form factors (valence quark contribution) which fits the data at intermediate and high scale, while 
the derived constraints at zero recoil [see Sec. II] are also satisfied. Then on top of the valence (or bare) form factors 
we place the meson cloud contribution which is relevant only at small and calculated consistently using effective 
chiral Lagrangian Finally we might conclude on the role of the meson cloud in the infrared domain. 

Due to the finite size of the source of the meson fields, meson loops are expected to be strongly suppressed for large 
Q^. To mimic the effect of additional regularisation of the integral, which leads to a restriction of the meson cloud 
contribution in the low region, we introduce the cutoff function fcut{t)- The modification is such that 

Fi™(t) ^ fcutit)F,1^l°-^{t). (97) 

The specific form of fcutit) should not modify F^2)°^'^{t) in the low region, but should diminish its contribution 

beyond a certain Q^. This function fcut{t) should have zero (or almost zero) slope at the point = GeV^ to 
ensure that it will not artificially contribute to the slope of the form factors. First the step function seems to be an 
appropriate choice for fcutit) but its sharp boundary affects the continuity of the form factors, and hence one must 
be careful. To avoid the problems associated with a sharp boundary we choose a smeared-out version for fcutit) with 

1 + expi- A/B) 

^™*(*^-l + cxp[(t-A)/B]' (^^^ 

where the parameter A characterizes the cut-off and B the smearing of the function. 

In our fitting procedure we use the experimental data (Refs. .58, 59, 60, 61, 62, il HI El El Elil El 113. El 
ElEllSlZSllZellTTllTllTliSililiiH ratios of the electromagnetic 

Sachs form factors of nucleons to the corresponding dipole form factors and on the neutron charge form factors. For 
the axial charge of the constituent quark g = 0.9 we first fit the contribution of the bare part, where the parameters 
are contained in the ansatz of Eq. H82() . Our parameter results for the fit of the bare part are (in units of GeV) 

A„i; = 2.0043, XdE = 0.9996, A„m = 7.3367, A^a/ - 2.2954, 

AuE = 0.8616, AdE = 0.9234, A„m = 0.9278, A^m = 1-0722. (99) 

Note that the parameters 7^ and 7^ entering in Eq. (|82|l cannot be fixed by considering only the bare part since they 
also occur in F(^ it) and F^ it) . The complete fit to the full data on the electromagnetic form factors of the 
nucleon fixes the remaining low-energy constants in the effective Lagrangian as 

7„ = 1.081, 7<i = 2.596, C2 = 2.502 GeV"\ C4 = 1.693 GeV"\ Jio = 1.110 GeV"^ 

67 = -0.650 GeV-^ eg = 0.030 GeV~^ ew = 0.039 GeV~^ , (100) 
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where in the current manuscript we put ky = for simpUcity. The coupUngs dio, 67, eg and eio are defined in Appendix 
C. With the choice A = (0.4 - 0.5) GeV^ and B = 0.025 GeV^ in f^^tit) the meson cloud contributions are suitably 
suppressed for values larger than t ^ 0.5 GeV^. Note, that from the analysis of the nucleon mass and meson- nucleon 
sigma-terms we get the following constraint on the parameter C2 = (1.093 — 1.354 §3 + 0.474 64 — 0.338 65) GeV~^. 
The best description of the electromagnetic properties of the nucleon is obtained at the value of C2 = 2.502 GcV~^ 
fixed from the data on the electromagnetic form factors. It means that we obtain the following constraint on the 
LECs 63, 64 and 65: 1.354^3 - 0.474^4 + 0.338 65 = -1.409. 

Results for the measured ratio of the electromagnetic form factors of the nucleon to the dipole form factor Gd {t) — 
(1 +t/0.71 GeV^)~^ are presented in Figs. 5-7. In Fig. 8 we present our result for the ratio of charge to magnetic form 
factor of the proton. For completeness, the electromagnetic form factors of the nucleon are also shown in Fig. 9, where 
the meson cloud contributions are shown explicitly in comparison to the dipole form factor. The role of fcut{t), which 
restricts the meson cloud contribution to the low region, is indicated in Fig. 10. Although the bare constituent 
quark contribution is fully parameterized, a consistent explanation of the form factors can only be achieved, when 
meson cloud corrections are included. We stress that the bare constituent quark Sachs form factors and hence the 
magnetic form factors at zero recoil {t = 0) are determined by the general constraints discussed in Section II. E. At 
this point the specific form of the quark form factors of Eq. H82|) . leading to a satisfactory description in particular 
beyond t ^ 0.5 GeV^, is not required. The results concerning the charge and magnetic radii of nucleons are 

= 0.881 fm, (r^)^ = -0.1177 fm^ 

= 0.869 fm, rlj = 0.847 fm. (101) 

The experimental values reported in Ref. 53] for the charge radii of nucleons are = 0.875 ± 0.007 fm and (r^)^ — 
-0.1161 ± 0.0022 fm^ For the magnetic radii, the analysis of Refs. [zl 113 gives r^,^ = 0.855 ± 0.035 fm and 
rl^ = 0.873 ± 0.011 fm, respectively. 

For further illustration we follow Ref. [H^l to deduce the radial dependence of the charge and magnetization densities 
of nucleons. The charge and magnetization densities of nucleons in the rest frame are 

p%{r)^- dkk^jo{kr)p%{k), pliir)^- dk jo{kr)pZ{k), (102) 

71" Jo Jo 

where fc^ — t/{l + tn) and jo{kr) is the Bessel function. The intrinsic form factors p^{k) and p^.j{k) are 

2N /n2\ 



p^(fc) = G^(g^)(l + r^)-, pj^,(fc) = ilMWi(i + ^^)-. (103) 

Pn 

Restricting to the discrete values — 0, 1, 2, where the full range is discussed in Ref. |5^, the results for the charge 
and magnetization densities of the nucleon are shown in Figs. 11 and 12. 

In Ref. | 55j the electromagnetic form factors of the nucleons are represented by the phenomenological ansatz 

G'^(Q') -G,(02) + ab.Q2G,(Q2)^ (104) 
where the "smooth" part Gs{Q^) and the structured or "bump" part Gb{Q^) are parameterized as 

^'^^ ^ (l + QV«ii)' (l + QVa2i)2' 

GbiQ^) = e n "6 J -i-e H "6 J . (105) 

The parameters aio, an, 020, 0,21, ab, Qb, and at are obtained by a fit to experimental data and the values are reported 
in Table 2 of Ref. 55]. The meson cloud part of our evaluation cannot be directly compared or matched to Gb(Q^), 
since meson corrections contribute to the magnetic form factors even at — 0, which is not the case in the treatment 
of Ref. [53. However, charge conservation restricts the meson cloud not to contribute to the charge form factors at zero 
recoil. Therefore, we can compare our results for the charge form factors to the phenomenological ones of Ref. |55j ]. 
Fig. 5 shows such a comparison for the case of the charge form factor of the proton. Based on the phenomenological 
ansatz of Ref. [H^ our result can be reproduced by readjusting the parameters of Eq. (|105|l as compiled in Table 3, 
which are not so much different from the original analysis of Ref. |55j . 
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D. Coupling of vector mesons to the nucleon 

Finally, we calculate the strong vector meson-nucleon form factors pNN and luNN at one loop. We follow the 
strategy already developed for the electromagnetic nucleon form factors. The corresponding bare operator is derived 
from the tree-level Lagrangian 

J^,yiQ) = J d'x e-^"- j^^^^x) , j^^^{x) = gv,, q(^^, + |^ a,. ^ 9 , (106) 

where qd'^ q = q{d'^ + d") q and Ay is the corresponding flavor matrix: Ap = diag{ 1,-1,0} for the /9° and A^^ = 
diag{l, 1, 0} for the w meson. The diagrams contributing to these quantities are displayed in Fig. 3: tree-level diagrams 
(Figs. 3(1) and 3(2)) and one-loop diagrams due to the dressing by a cloud of pseudoscalar mesons (Figs. 3(3) and 3(4)). 
The Fourier transform of the dressed vector-meson-quark transition operator has the following form 



4T(i) = I I d'^e-^^^qix) 



q{x) (107) 



where fv^ and fvp are the matrices of the Dirac and Pauli form factors describing the coupling of u, d and s quarks 
to the p and lo vector mesons. These form factors are given in terms of the Euclidean values of momentum squared 
t — —q^ with: 

fVoit) ^ J2 fvnit)^ fvAt)^ E /Vp(0, (108) 

i—u,d.s i—u,d.s 



/;p(i) = -,C(i) = kv 9v,, { 1 - r<o(i) + :j"h'o(^) \ , 



1 
3 

/:V(0 = fUt) = kvgv,, |l + 3rr4o(t) + im'^(t)| , 

fpp (t) = 0, f^^ (t) = kv 9vn I -4 mfo (0 } , 

where A^ = diag{0,0, 1}. Here (t) and m^{t) are the meson-cloud contributions given in Appendix C. 

To project the dressed quark operator p07ll onto the nucleon we proceed in analogy to the electromagnetic operator 

( N{p')\ J^:ril) \N{p)) = {2-Kf 5\p' ^p^q) u^ip') 1 (7^ GvNNiq^) + ^ a^.q" FvNN{q')]uNip) 

= {2nf 5\p' -p-q)\ |/{/Jg2) (^(p')l j':S'(0) \N{p)) + i ^ flil^) {N{p')\jll%{Q) |7V(p))| , (109) 

where the bare matrix elements (.^V(p')| j^jJ™-(0) \N[p)) are defined in Eq. Here Gvnn{<1^) and FvNN{<f) are 

the vectorial and tensorial couplings of vector mesons to nucleons. We can express the strong pNN and loNN form 
factors through the bare electromagnetic nucleon form factors. After a simple algebra we arrive at: 

GpNN {t) = 9v,, [ Ff- (t) (t) + Cr (t) Fp, {t) ] , 

G^NNit) = 5v,,[Fr(t)L",(o + Gr(t)/:,(t)], 

FpNN it) = gv,, [ Fr (t) (t) + Gr it) it) ] , 

F^NNit) = gv,, [Ft it) f:jt) + Gl+it) fZ.it) ] , (110) 
where t ^ -q^ and = i/f" ± uf with H = F ot G and / = 1 or 2. 
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Finally, we present the expressions for the values of the vector-meson nucleon form factors at zero recoil or the 
coupling constants Gvnn and Fvnn which originate from the nucleon- level Lagrangian (for details on the nucleon- 
level Lagrangian see Ref. 

CvNN = \N\(r GpNN - o^" p^T+ ("7^ G^NN " a^'' d.) cj An. (Ill) 

After a simple algebra we arrive at: 

/~i „ PpN N bare ,,barc 1 , i / rbarc rbarc \ ri 1 r 1 

GpNN = gvqq , = ^J■p - M« - 1 + fcy ( d - ) [I + dpi , 

UpNN TT^q 

G^NN = 3 gv,, , = Mp^'" + - 1 + fcy — ( + ) [1 + S^] , (112) 

^ulNN TTlq 

Sp = -m^„(0) + im?o(0) , = 3m?o(0) + im^„(0) (113) 

are the corresponding one-loop corrections. For fcy = these equations reduce to the well-known SU(3) relations, 
relating the matrix elements of the vector current with different flavor content ^9?\ . Using the actual expressions for 
the bare magnetic moments and the tensor charges of the nucleons we finally obtain for the ratios of the tensor and 
vector couplings: 

FpNN _ f 9a\^ f -. , n I r A 1 F^NN _ 1 f QA^ ^ 



where 



„ , , ,l + fcy[l + <5,] -1, ^ 1 + 3fcy[l + ^c.] -1. (114) 

^pNN \ 9 / \ J ^ujNN i) \ g / \ ) 

Note, that numerically the one-loop corrections to the tree-level results for the strong vector-meson nucleon form 
factors FpNN and F^nn are rather small: 5p = 0.005 and S^j — 0.011. 



V. SUMMARY 



We developed a manifestly Lorentz covariant chiral quark model for the study of baryons as bound states of 
constituent quarks. The approach is based on the effective chiral Lagrangian involving constituent quarks and the 
chiral fields as effective degrees of freedom. This Lagrangian is used in the calculation of the dressed transition 
operators which are relevant for the interaction of quarks with external fields in the presence of a virtual meson cloud. 
Then the dressed operators are used in the calculation of baryon matrix elements. 

Our main result is as follows: we perform a model-independent factorization of the effects of hadronization and 
confinement contained in the matrix elements of the bare quark operators and the effects dictated by chiral symmetry 
which are encoded in the corresponding relativistic form factors [see e.g. Eq. (I19|l]. Due to this factorization the 
calculation of chiral effects and the effects of hadronization and confinement can be done independently. All low- 
energy theorems are reproduced in our approach due to the chiral invariance of the effective Lagrangian. The evaluated 
meson-cloud corrections are in agreement with the infrared-singular structure of the corresponding nucleon matrix 
elements In particular, we reproduce the leading nonanalytic (LNA) contributions to the nucleon mass, to the 

pion-nucleon sigma-term, to the magnetic moments and to the charge radii of the nucleons. The LNA contributions 
to the nucleon mass and magnetic moment are proportional to the Mp and Mp, respectively, where Mp is the 
pseudoscalar meson mass. The nucleon radii are divergent in the chiral limit. Using model-independent constraints 
on the bare constituent quark distributions in the octet baryons, we work out model predictions for the magnetic 
moments. Based on a full parameterization of the bare constituent quark distributions in the nucleon, we give results 
for the full momentum dependence of the electromagnetic form factors of the nucleon and indicate the role of the 
meson cloud contributions. Presently, the calculation of the matrix elements of the bare quark operators is, besides the 
model independent constraints, based on parameterizations. The direct calculation of these matrix elements should 
be performed in quark models based on specific assumptions about hadronization and confinement. 
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APPENDIX A: CALCULATIONAL TECHNIQUE OF LOOP DIAGRAMS 



The calculational technique of the loop diagrams in Figs.L 2, 3 and 4 is referred to as the infrared dimensional 
regularization (IDR) and has been discussed in detail in Refs. |5lllCl||. Our case here differs by the use of the constituent 
quark degrees of freedom instead of the nucleon. We briefly recall the basic ideas of this technique. As we discussed 
before, in Baryon ChPT loop integrals are non-homogeneous functions of the mesonic momenta and the quark masses 
due to the presence of a new scale parameter, the nucleon mass. As a result, the loop integral contains an infrared 
singular part containing the fractional powers of the meson masses and an infrared regular part involving the fractional 
powers of the nucleon mass. Due to the presence of the regular parts in the loop integrals, their chiral expansion starts 
at the same order as the tree graphs, it therefore spoils the power counting rules. The idea of the IDR method is to 
remove the infrared regular parts of the loop integrals from the consideration and to absorb them in the low-energy 
couplings of the underlying chiral Lagrangian. The IDR method is consistent with Lorentz and gauge invariance. 
Also, the chiral power countings are preserved and the Ward identities of chiral symmetry are fulfilled. 

The method can be explained by using the simple example of the self-energy graph shown in Fig. 4(3) contributing 
to the quark mass operator We consider the scalar loop integral in d dimensions: 

^^P^^ J (^;;^ [M2 - fc2 - ie] [m2 - (p - fc)2 - ie] ^^^^ 
where M and m are the meson and constituent quark masses, respectively. Using the master formula in d-dimensions 

d'k _ r(n + d/2)r(m-n-d/2) ,(„_„)^^ 

(27r)'i [M2 - A:2]™ [At^Y/^ T{d/2)V{m) ^ ' 

we get at threshold ^ {m + M)"^: 

rr, 2^ M<^~^ m'^-^ r(2-d/2) 

m + ili m + Ai (47r)'^/^ (a — 3) 

=/(p2) =fl(p2) 

Here /(p^) is the infrared singular piece which is characterized by fractional powers of M and generated by the loop 
momenta of order of the meson mass. For /(p^) the usual power counting applies. Another piece in the decomposition 
of H{p^) defined by R{p^) is the infrared regular part which is generated by the loop momenta of the order of the 
constituent quark mass m (in our counting m is of the order of ^xSB ^ 1 GeV). As discussed before, we remove the 
regular part R{p'^) by redefining the low-energy coupling constants in the chiral Lagrangian. In Ref. ^ a recipe was 
suggested how to split the integral H{p^) into the singular and regular parts. We use the Feynman parametrization 
to combine two multipliers a = Af2 ^ k"^ — ie and b — rn^ — (p — fc)2 — ie in the denominator of H{p^): 

d'^k 1 r d'^k } dx , 

(A4) 



{2TTy ab J {2'kY J [a{l-x) + bxY 



and then to write down the integral from to 1 as the difference of two integrals: 



dx 



dx ... (A5) 



Then the integral from to oo is exactly the infrared singular part and the integral from 1 to oo is the infrared 
regular one. This method can be applied to any general one-loop integral with adjustable number of meson and quark 
propagators. The calculational technique suggests that: one should separate numerator and denominator; simplify 
the numerator using the standard (invariant integration) methods. Finally, the result can be reduced to the master 
integral I{p^) and its derivatives (like in the conventional dimensional regularization). The integrals containing only 
the quark propagators do not contribute to the infrared singular parts, and therefore, vanish in the infrared dimensional 
regularization. 

The ultraviolet divergencies contained in the one-loop integrals are removed via the renormalization of the low- 
energy constants in the chiral Lagrangian. To perform the renormalization of the constituent quark at one loop 
and to guarantee charge conservation we need the Z-factor (the wave-function renormalization constant), which is 
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determined by the derivative of the quark mass operator (generated by the diagrams in Fig. 4) with respect to its 
momentum: 



9So 



q = u,d,s . 



The matrix Zq = diag{Z„, Z^, Zs} with Z^ — Zd = Z is given up to fourth order by 



Zq = I 



P=n,K,ri P 



g ap Ap + Qp 



where 



2A/|, 



Qi 



247r2 



A(/i) 



ap 



167r2 \ ^l ) 



A(m) 



(i-4 



-1 + 



37r Mp 



3 Mj, 
2 m? 



(47r)2 
647r2r7i 



1 



OLK 



-3Q 



2/+ -A3 , 
2 ^ 



-4q 



1 3, 

2 + 4^- 



(A6) 



(A7) 



-(ln47r + r'(l) + 1) 



(A8) 



4 1 
/3,r = 1 , = 3 ' ^ 3 ■ 

For the evaluation of the form factor ff){q^), the quark charge Q [diagram in Fig. 1(1)] has to be multipHed by Zq. 
For fp{q^), a second-order contribution to the quark anomalous magnetic moment proportional to the cg [diagrams in 
Figs. 1(2) and 2(2*)] has to be renormahzed by the Z^-factor. The term proportional to C2 must be dropped, because 
the product cg C2 is of higher-order when compared to the accuracy we are working in. 



APPENDIX B: LOOP INTEGRALS 



In this appendix, we present the loop integrals contributing to the electromagnetic transition operator between 
constituent quarks and evaluate them in the infrared regularization scheme i. Originally, these integrals have been 
introduced in chiral perturbation theory (ChPT) [H-flol]. 

For the momenta of initial, final quark and photon field we introduce the notations: p, p' and q = p' —p, respectively. 
We also introduce the notation P = p' +p. Since the external quarks are on the mass shell p^ — p'^ = m? , the structure 
integrals can be expanded through a set of scalar functions which exclusively depend on the transverse momentum 
squared t ^ — —q^ and the masses of meson M and constituent quark m. As usual 0,0], we identify the masses 

o 

of particles occurring in the loop integrals with their leading order values, Mp —^Mp and ruq — > m. For universality, 
we calculate the integrals for adjustable values of the constituent quark mass inside the loop (m*) and for the external 
line (to). Finally, in numerical calculations we will neglect the difference between the masses m* and m, setting 

TO* = TO. 

1. Infrared parts of loop integrals 



We deal with the following loop integrals: 
d'^k 1 



(27r)'^i AP ~ fc2 
d'^k 1 



{21:)^ [M2-j!c2][jVf2_(fc + g)2 



(Bl) 

j(°)(i,M2), (B2) 



d'^k k 



-^2 n , ^21 = -lq,j'^'\t,M'), (B3) 



{2TTYt [M2 - fc2] [M2 _ (fc + g)2] 2 
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[M2 


-A;2][M2- 


-{k- 




d<^k 




1 








[Af2 


— /c2][to*2 


~iP 


-fc)2] 
















— A:2][j7j*2 


-iP 


-m 


d<^k 








1 




[Af2 


— /c2][^*2 




-A:)2][ 


d'^k 











AO) 



= I'^^>it,M',m',m*') 



= P /"(^Vf M2 ?7i2 7J7*2\ 

(27r)'ii [M2-fc2][TO*2_(p_^.)2][^,2_(p, _fc)2] ^ 12 l^, , m , m 



g^,/^^)(^,M2,m^m*2) 



f Jl. 

J (2^ 



{2'KYi [M2 - )!c2][to*2 _ _ /j)2] [m* 2 _ (p/ _ fc)2j 

+ P^P,/(^)(t,M2,m2,m*2) + Q^g,/(^)(t,M2,m2,m*2), 

i = m2 m2 m*2^ 

(2^)'i2 [M2 - A:2] [Af 2 _ (fc + ^)2] 2 _ (p _ ^)2] 21 I > > " : J . 



(27r)'^i [M2 ~ /c2][M2 - (fc + g)2][m*2 - (p - kf 



1 



q^4'i\t,M-',m',m*-') 



2 

d<^k 

{2-KYi [M2 - /c2][Af2 _ (fc + g)2][rri*2 _ (p _ /j)2 



^/J, ky 



(B4) 
(B5) 
(B6) 
(B7) 
(B8) 



(B9) 
(BIO) 



(Bll) 



(B12) 



where the symbol / represents the loop integration according to the infrared dimensional regularization scheme |5Llin|| 

2. Reduction formulas for loop integrals 

Higher-order tensorial integrals can be reduced to the basic scalar integrals using the invariant integration ^3 '■ 

1 



A{d~ l)t 



(t + 4A/2)j(°)(t,M2) - 2 Am 



/W(M2,m2,m*2) 



/«(t,M2,m2,m*2) 



4 ' ' ' 2t 
1 



M ■ 



2m? 



1 



4rn2 -|_ t 



M*^I<^^\M^,m^m*^) + Am 



(B13) 
(B14) 
(B15) 
(B16) 



d-2 



1 



(d - 2) (4m2 + t) 



M*^d-l)l[l\t,M\m^,m*^) 



1 



{d-2)t 



-M*^l\l^{t,M\m\m*^) 



(2M*2 + t)4f{t,M'\m\m*'') 



2 (4to2 + 1) 

(4M2 + ^)4;^(^,M^m^m*2) 



4(d- 2) 



2(2M*2 + t)4l>{t,M\m^,m*^) ~ 2 /^(M^, m^, m* 2) 



-(4M2 + t)4?(t,M2,m2,m*2) 



4 (d - 2) (4m2 + 



2(d-l)(2M*2 + f)4J)(t,M2,m2,m*2) 



+ 2/(")(M2,m2,m*2) - 2 (d - 2) /(^^ (M^, m^, m* 2) 



-(4M2 + (d- l)i)4°^(t,M2,m2,m*2) 



4(d- 2) 
2(2M*2 

2(d-3)/(")(A/2,m2,m*2) + 2 (d - 2) /(^^ (M^, m^, m* 2) 



2(2M*2 + t)4J^(t,M^m^m*2) 



where M* ^ ^ M'^ + m'^ - m* . 



26 



(B17) 



(B18) 



(B19) 



(B20) 



(B21) 



(B22) 



(B23) 



3. Scalar loop integrals 



The scalar loop integrals are given as [slITol 



Am = 2M^\m, 
J^°\t,M^)^-2\M 



l[l\t,M^,m^,m*^) 



Am — 



M 



d-4 



(47r)2 {d-A 2 



^ i(ln47r + r'(l) + l) 



1 



167r2 

M*2 



A A/ — 



t 



32 7r2 



32 7r2 



32 7r2Mm3 "^177^2 



g( -!^,M,m,m' 



(B24) 
(B25) 
(B26) 

(B27) 
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Am + 



32 7r2 



+ 



1 



32 7r2 m2 



/ill j^,M,m,m* 



■2l jj + ^) /i2(^,M,m,m^ 



where 



2m*M ^ m 



M 



M = 



r = 



m 
m 



The dimensionless functions k, f, g, hi and /12 are given by the expressions 



fc(s)= / dx\n{l + x{l-x)s) = J^ in^Si^-2, 



i 

f(s) = [ dx ^ 

J l + x{l 



In 



x)s y's(4 + s) V4 + s - ' 

i/r* 2 + a:;(l — x),s 



arccos 



g{s, M,m,m*) = / c?x 



/ii(s, M, m, m' 



-r* (f^ + M*)- 



1 + a;(l — a;)s 



[1 + - x)s)] ^{l-n'^)r*^ +x{l- x)i 



)= / dxj 



ln[l + - x)s] 



+ jj, x{\ — x)s 



1 

/12 (s, M, rn, m*) = J dx arccos 



[li*{l + x{l-x)s) +n] r* 



\/\ + /Lt^ — x) s ^1 + sa;(l — a;) 
1 



[1 + ii^xiX- x)s\^J\ + 02 + a;(l-a;)s 
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We exphcitly show the contributions of the diagrams in Figs. 1 and 2 to the Dirac and Pauli quark form factors. 
We use the following notations flj\t) and fp\t) where the superscript i refers to the numbering of the diagrams. 
We expand the diagonal 3x3 matrices fjj{t) and fp\t) in the basis of the 3x3 matrices Q, I, and A3: 



2/3 0\ /100\ /lOO 

Q=\ -1/3 ,/= 010 ,A3- 0-10|. 

-1/3 J yooiy \ooo 

The contributions of the diagrams in Fig.l to the Dirac form factors of the u, d and s quark are 

fS\t) = fo^it) = fS\t) = fl^'\t) = 0, 
f^^\t) = 2tdwQ, 



(CI) 



J2 Afef, 1 = 5,6,7,9,10,11, 

P=w,K,ri 



(C2) 



where Zq is the finite part of the renomahzation matrix Zq defined in Eq. ljA7|l . 



A? 

^5 


- O + 




---Q 
" 3^ 


9 


Afi 


= A? = 

7 




- Af = 

7 


= Af 

9 






+ i/-A3) 


^10 - 










^11 — 


3^ 



i/-iA3. 
9 3 ^ 



3^9 3 



9 3 



= -- 



4 



Ff 



1^ {t, Ml) - Aw?I^^^ {t, Ml, m^^m^) - lQm^Pi> {t, M^, m^,m'^) \ , 



^ P 



-'j(^){t,Ml,m\m^), 

2tg^m'^ -(3) 2 2 2n 
-^^Il2 {t,Mp,m ,m ), 

3C2 



/V/r4 



The contributions of the diagrams in Fig.l to the Pauh form factors of the u, d and s quark are 



/p\t) 



p(7) 



(8). 



-2tdwQ , 
—AmteioQ — 16m 



Af , 

P=Tr,K,q 



M^eeQ + -{Ml - Ml){e7[Q - A3 - //3] - eg/) 
i = 5,6,10,11,12, 



where 



A12 — A3 , A^ _ 3(5 ~ A3 , A^2 — , 



^-^|/(i'(M|>,TO2,m2) - Mlh2{t,Ml,m^,m^)^Ap^^{t,Ml,n?,'m^ 



Fl 



+ 2t (^7{2^(t,M|,,m^TO2) -/j2^(t,M|,,TO2,m2) 
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mfi 



p 



(C7) 



The contributions of the diagrams in Fig. 2 (vector mesons contributions) to the Dirac and PauU form factors of the 
u, d and s quark are 



fV{t) = f^''\t) = f^^'\t) =0, 



t 
2 

kv 



+ 



9 

kv 

fp '{t) = fp'it) =0 



ip 



Z Xp Dp{t) + ^ZX^D^{t)+'^ ZsX^D^it) 



+ ^[3XpDp{t) + X^DUt) + 8X^D^{t)]m1o\ , 



f^P^'\t) = -lkv\XpDpit) 



p2' 



' ji^)it,Ml) + ^ji^){t,M%) 
+ ^D^{t)J^'\t,Ml) + ^I?,(t)j(i)(t,Mi)l , 



where Ap = As, = 2Q + 2//3 - A3 , A^ = 2Q - 7/3 - A3 and = 1/(M^ + i). 

In the expressions above, the convergent parts of the structure integrals are 

1 



l[l\t,M\m\m^) = 
l[l\t,M^m^m^) = 



5767r2t 



[t + 3k{t) {t + 4M2)] , 



647r2 
1 



327r2m2 

,3 



2567r2 
D 



[H - g{0,M,m,m){4-fi^)], 
[g{t,M,m,m)iJ,- f{t)], 

-2ij,Dr + g{t,M,m,m){'i- ii^Dr) + g{Q,M,m,m){A- j/)Dr 
2/i(l + 2£)^) - iJLf{t) - g{t,M,m,m){4:-3n'^Dr) 



10247r2m2 ^ 

5(0, M, m, m) (4 - i^^){2 + SD^) 



l[^{t,M\m\m'') = 



UrDr - 2g{t, M, m, m)(l - /i^Dr 



I\l>{t,M\m\m') = 



- g{0,M,m,m){8- 4:1^^ + {12^ - 4)Dr) 
1 



2Dr/i{n' - 2)(2/x - g{0,M,m,m){A - fj,')) + k{t){2 + {/i' - 2)D^) 



I^;>{t,Ml,m\m'') 



5127r2 ^ 

+ 2Hi2{t,M,m,m){4- {2- n^fDr) 

Dr 



10247r2m2 



2(4 + - + 2(m^ - 2)(m^ + 1)-Dr) - 6fc(i)(2 + (/x^ - 2)£>^) 



(C8) 
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+ 2/(t)(2 + - 2)Drf + 5(0, M, m, m)M(4 - m')[2(m' + 2) + 3{ii^ - 2)Dr] 



Hi2{t, M, m, m) (12 - S/x^ + 2r - - 2fDr) 



(C9) 



where 



Hi2{t,M,m,m) = hi{t,M,m,rn) + h2{t,M,m,m) 
1 



(CIO) 



1 + t/4 



In the practical calculation, we keep the exact form of the function f{t), while we expand the functions k{t), 
g{t, M, m, m), hi{t, M, m, m) and /i2(i, M, m, m) in powers of 9 and ^ with 6 = t/M"^ and = M/m. 
Finally, we keep our expressions for ef , and mf up to 0(^p, /ip) which are explicitly given as 



4 



327r2F2 L 



/(t) 1 - ^M^^D^ + ;,^p 1 _ _ 2Z?^ + ^ (1 - - r i^. 



/4, 



2 2 



A^p ( 1 - - ) + T 



2 

7 357r 13 2 1057r 3 107 4 



+ ^Mp 



(Cll) 



£7 = 



g'^m'^lip 
167r2i^2 



1927r2F| ' 



^10 



7r(l + 2D,) + ^ I 1 + ) + -apDrfit) 



^11 



3C2TO Vp 
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g'^m'^li^p ( 3 



167r2F2 I 8 



2 2 r 

g m 



^lipf{t)Dl + D. 
lip 



3/2 ^ 
8 I3 ' 



15 2 4 3 

yTT/ip - -lip 



+ 



96/ip 



TT + Slip 



IOStt 



-Aip - yA*p 



52 3 15757r 4 



107 



12^^^^ + To"^^ 



m 



10 



222 
g m Hp 

327r2F2 



/(O 



l--i^,A'p(l-^r) 



TT/XP ( 1 + -7^, - 2i^^ 



2 ^ 



m 



11 



3c2m lip 



p CiUlt 

After absorbing divergences coming from vector meson-exchange diagrams the renormalized LECs are written as 
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" = '^'-ffi^ff '"-^• 
^. -»M-5^inf (C12) 



where the /3-coef5cients with taking into account of vector mesons are given by 

1 + 552 



4 



and 



/3eio = 7r~ (1 ^ '^5^ + kv + 417104) , 
8m 

= -77;;— C6.g + kv) - -T— C6 - 4C4TO - — , 

144m lorn 9m 

/^e7 = ~7^(c6 + fcv) ff^ - (c6 - 4c4m) - , 
48m lom Dm 

Pes = -T^{c6 + kv)9^ + -^ic6-4.Cim) + -f- (C13) 
lom lom 4m 

APPENDIX D: MESON-NUCLEON SIGMA-TERMS 

Two important identities are involved in the evahiation of the meson-nucleon sigma-terms. In particular, the 
derivatives with respect to the current quark masses are equivalent to the ones with respect to the meson masses |4^ : 

d ^ , /I 9 2d 



- ^ = AfiJ-^ + ^-^lS,, (D4) 



with = -M|. + 3M2/2 . 

Below we give the exact expression for the typical type of derivatives: 



9m 647r \ 3 "^7 F^^ j 327T^m{F^ 3F^ 27 F^ 

+ 2ei A//2M2 - (Af^ - A/^) + ^e^M^ " g Af^ _ ^^2) ^ 



One should note that there are additional useful relations between different sigma-terms 22]. In particular, with the 
'kn 



definitions of j/jy E^nd crj^^ we can relate KN and rjN sigma-terms to the ttN sigma-term as 



m + m^ 7=1 d u ryT^i „ „ 7h + 2yNms 



4m 6m 
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Table 1. SU(6) couplings a^* and a 
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Table 2. Magnetic moments of the baryon octet (in units of the nucleon magneton /ijv) 
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Table 3. Comparison of the parameters of the original phenomenological ansatz in Ref. to an equivalent set, 
which is obtained by fitting to our results for the charge form factor of the proton. 
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Fig. 1. Diagrams including pseudoscalar meson contributions to the EM quark transition operator up to fourth 
order. Solid, dashed and wiggly lines refer to quarks, pseudoscalar mesons and the electromagnetic field, respectively. 
Vertices denoted by a black filled circle, box and diamond correspond to insertions from the second, third and fourth 
order chiral Lagrangian. 



37 



(2*) 



V 



(3*) 




Fig. 2. Diagrams including vector-meson contributions to the EM quark transition operator. Double-dashed lines 
correspond to vector mesons. The symbols V and T refer to the vectorial and tensorial couplings of vector mesons to 
quarks. 
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Fig. 3. Diagrams contributing to the quark operator describing stong interaction of vector mesons to quarks. The 
symbols V and T refer to the vectorial and tensorial couplings of vector mesons to quarks. 
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Fig. 4. Diagrams contributing to the mass operator of the quark at one loop. Vertices denoted by a black filled circle 
and diamond correspond to insertions from the second and fourth order chiral Lagrangian. 
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Fig. 5. Ratio G^^{Q^)/Gd{Q^) '■ (o-) Overall range, (b) Up to — 2 GeV^, the solid line is the total contribution, 
the dotted line is the bare contribution and the dashed line is the result reported in Ref. Experimental data are 

taken from Refs. \M,M,\M,M,\M,\M, EI, M]- 
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Fig. 6. Ratio G\i{Q'^) / {upGoiQ'^)) • (a) Overall range, (b) Up to Q'^ = 2 GeV^, the solid line is the total 
contribution and the dotted line is the bare contribution. Experimental data are taken from Refs. \6(]i . \ 6A \ 6fl. \6'i. f^. 
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Fig. 7. Ratio G'^'^^(Q^)/(^„ G_d((3^)) ; (a) Overall range, (h) Up to — 2 GeV^, the solid line is the total 
contribution and the dotted line is the bare contribution. Experimental data are taken from Refs. f?^ . It^. ItI. It^. [tJ 
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Fig. 8. Ratio UpG^^iQ"^) / G\.[{Q'^) in comparison to the experimental data. Experimental data are taken from 
Refs. ^,[11 
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Fig. 9. The charge and magnetic form factors of the nucleon and the contribution due to the meson cloud. The 
corresponding dipole fit f or G^tp ( Q^), Gj i AQ '^) and G^{Q'^) are also presented. Experimental data for G'%{Q^) are 
taken from Refs. ^\M,\M,\MA82,\M,\M,\M, I12/- 
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Fig. 10. The meson cloud contribution to the nucleon charge and magnetic form factors with the cutoff Junction 
fcutiQ'^) (solid line) and without fcut{Q^) (dotted line). 
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